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Abstract 

Landau's theory of electron motion in stationary magnetic fields is extended 
to the inclusion of bouncing along the field between mirror points in an 
inhomogeneous field. The problem can be treated perturbation theoretically. 
As expected, bouncing is quantized, causes a weak shift in Landau levels, 
reduces parallel degeneracy, but does not contribute to diamagnetism. 
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In a homogeneous magnetic field B = Bqz, with vector potential A(x) = 
-Bqu, 0, 0) in the Landau gauge, electrons perform Landau-cyclotron orbits 



turn mechanical problem had been treated long 



Uc = eBo/nie. Tl 
ago by Landau ( 


le quan- 


Landau 


Landau & Lifschitz 




1965 



Kleinert 



2004 . and others) and yielding the celebrated quantization of per- 



pendicular electron energy e± = fvjJc{L + i) into the chain of Landau lev- 



els L £ N an d theory of electronic diamagnetism (cf., e.g.. 



Landau 



1980 : 



Huang . Il987l ). Clearly, in homogeneous magnetic fields, the parallel energy 
eii = p1/2mf, = h^k1l2me is unaffected. Since the problem is independent 
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on y,Py, Landau levels are highly degenerate. In a homogeneous field their 
degeneracy is g = r^Vz/X'^^, where rx,rz are the transverse to y extensions 
of the spatial volume under consideration, i.e. the perp endicular surface ex- 
press ed in terms of a magnetic length = ^/27rh/eB() ( Aharonov Sz Bohnj . 



19591). 



In a magnetic mirror geometry (like, e.g., the Earth's radiation belts) elec 
trons of sufficiently low parallel energy may become tr apped, bouncing be 



tween the mirror points along the r nagne tic field (cf., e.g.. lKivelson fc Russelll . 



1995 



Baumjohann fc Treumannl . 



19961 ). Classically this motion is under- 



stoo d as an oscillatio n along the magnetic field at bounce frequency Ub (cf. 



e.g., 



Hasegawa . 



19751 ) . This frequency is formally given by the inverse double- 



bounce time Th which, for an electron of energy e, is an integral over the 
parallel energy e\\{s) = e — e_L(s) along the magnetic field: 



n = 



27T 



2mp 



S2 



ds 



Vi-e±(3)A 



(1) 



between the two mirror points si,S2- These are defined as the points along 
the magnetic field, where all the electron energy is in the perpendicular, 
i.e. e||(si^2) = 0. For non-relativistic electrons, neglecting radiative losses, 
the electron conserves its magnetic moment /i = e±{s)/ B{s), and we have 
e±(si,2) = e, or e±{s)/e = B{s) / B{si^2)- 

If we assume that the magnetic field possesses some approximate sym- 
metry with respect to s (for instance like a dipole field) and that the mirror 
points are located not too far away from the symmetry plane, then si = 
—S2 = —Sjn, and -B(O) = Bq is the minimum of the magnetic field along the 



field. Expanding B{s) around minimum, B{z) ^ Boll + az"^ 



^B"(s) 



\s=0 



is the second derivative along B taken at field minimum. Close to the plane 
of symmetry s ^ z becomes a straight coordinate —Zm<z<Zm that is lim- 
ited by Zm = \/ {R — and R = B{sm)/BQ is the magnetic-field mirror 
ratio, the latter assumed being a known quantity. 

With these approximations the vector potential of the mirror field be- 
comes, in the Landau gauge, 

A[x, z) = [-BQ{l + az'^)x,0,0], \z\<Zm. (2) 

The correction introduced by the mirror geometry is small, and the vector 
potential still possesses only one component Ay which, however, now depends 
on the two spatial coordinates x, z. This complicates the problem and, in 
addition, gives rise to a weak transverse magnetic component By — dzA^ — 
—2Boayz. It causes a weak particle drift in y-direction which, in higher 
order, would ultimately break the Landau degeneracy. 

With the above vector potential the Hamiltonian of the gyrating and 
mirroring electron becomes 

+ pI+pI]- (3) 

Since, near field minimum Bq, a and z can be assumed small numbers, this 
expression can be substantially simplified by expanding the second term in 
the braces, yielding 

n^nL + v{y,z). (4) 

Here Hl = (2me)"^|p^-Fp^+"^e'^c(2/-2/o)^|, with yo = p^/nieUJc, is the usual 
harmonic oscillator Landau-Hamiltonian which (to lowest approximation) 
depends only on spatial coordinate y, and 

V{y, z) = -Ucaz^y (p^ - rUeUcy^ (5) 
3 



2m„ 



P:c - e-Bo(l az^)y 



is a perturbation which is considered to be small, a well justified assumption 
near field minimum along the field line. In the last two expessions px, y, z 
are operators. Strictly speaking, Uc{z) = cBq (1 + az"^) /me in T-Ll is a z- 
dependent cyclotron frequency. However, for the assumed small a, z values, 
it varies only very weakly along the field. Including this 2;-dependence would 
account for a slight shift of the Landau levels and changes in the Landau wave 
functions along the field. In the following it will, for our limited purposes of 
investigating the more interesting effects of bounce motion, be taken as the 
constant value Wc(0) = eBo/nie in the field minimum. 

Under these assumptions, the system is, by definition, periodic and thus 
stationary. With the Landau-Hamilton operator in the Schrodinger equation 

: bed Landau solutio n \L(x ) ) = \L ^ ^^x)) 



is known (see, e.g.. 



)), the unpertur 


Landau. 


1930 


1 



Landau fc Lifschitzl . Il965l : iKittel , 



1963 . 



and others), and the whole problem can be treated by perturbation theory. 
We note that the perturbed Hamiltonian remains to be independent of x. 
Hence the corresponding momentum may be replaced by px = hkx, and 
degeneracy in x remains unresolved even for a bouncing particle. Solving 
the perturbation problem will provide the energy shift = — e^"'' in the 
Landau levels; it also provides the perturbed wave function \L^^'^) ^ l-^^*^-*)- 

Before proceeding, we instead determine the corresponding quantities for 
the pure bounce motion. This can be done by realizing that the Hamiltonian 
Eq. ([3]) can be averaged over the fast Landau-cyclotron oscillations. This 
is easily done when writing classically for the two coordinates {x(t),y(t)} = 
rc {cos Wet, sin Wet} with the classical gyroradius. Correspondingly, the 
momenta become {px^Vy\ = meWcrc {— sin Wct, cos Wct}. Inserting into Eq. 



4 



([2D and averaging with respect to time over two gyroperiods Tc = 2n /uc 
(because the Landau ground state frequency is ^Uc) yields for the gyro- 
averaged pure bounce-Hamiltonian up to second order in z 

2 

nt = ^+^Bo{l + az^), (6) 

which is the Hamiltonian of a harmonic oscillator of reduced Hamiltonian 
= Tib — /i-Bo and energy eigenvalue e' = e — /i-Bq, with fi{L) = e^{L)/ Bq 
the magnetic moment of the electron in Landau level L. Its Schrodinger 
equation is "Hyl^) = {(i + \)tujj\\\i), which identifies 

4(L)= + (7) 
as the energy spectrum of the bounce oscillation in Landau state L, and 



-II W = . = [L + 1) « c.. (8) 

y me \ me 

as the corresponding quantum mechanical bounce frequency. The waye func - 



tion \i) can be found in any QM textbook (see, for instance, ISakurai Il994j ). 
As expected, the parallel energy of an electron in any Landau state L bounc- 
ing along the magnetic field is quantized. The restriction on the parallel 
quantum number is i/L < uJc/uj\\ imposing an upper bound on £ for every 
Landau quantum number L. Since Uc ^ uj\\, any perpendicular Landau level 
contains many parallel bounce levels. This resolves the parallel degeneracy 
reducing Landau degeneracy to g = r^/ Xm- Though this is not a particularly 
deep insight, it clarifies the properties of electron motion in mirror magnetic 
field geometries. 

One may also ask whether the bounce motion in an electron gas of density 
N introduces any contribution Xb to the diamagnetic susceptibility x- This 



correction is obtained from the second derivative of the logarithm of the 
bounce partition function logZh with respect to the magnetic field B at 
constant temperature T, volume V, density A^, a nd chemical p otential /2. It 



can be shown by conventional methods (cf., e.g. 

log Zb ^ A^V" 



Huang 



19871) that 



1 



1 / hoj 

24 



T 



(9) 



Inserting for one immediately finds (d'^logZb/dBQ) \t,v,p. = 0, hence yield- 
ing = 0. To lowest order the bounce motion does not contribute to mag- 
netic susceptibility. This is reasonable, because contrary to Landau orbits 
bouncing electrons do not cause a net current, whose magnetic effect thus 
vanishes. 

Returning to our initial problem, we may now determine the perturbation 
of any of the Landau levels caused by the parallel electron motion when the 
electron is forced to perform a bounce along a magnetic mirror field geometry. 
We are mainly interested in the energy shift 



f{L) 



(10) 



of a particular Landau level caused by the bounce motion. This can, to lowest 
perturbation order, be calculated si mply from the known unp e rturbe d Lan 
dau wave function see, e.g.. 



Landau fc Lifschitz 



1965 



Kittel 



Kleinert 



1963; 



2004 the latter for its path-integral amplitude) using the pertur- 



bation Eq. (|5]) of the Hamiltonian. The Landau wave function consists of 
products of exponentials and Hermite polynomials. Exploiting the orthonor- 
mality conditions in ?/, the term proportional to ypx in ^(z/? z) is linear in 
y\ being odd it vanishes by y-integration. The contribution comes from the 



6 



quadratic in z and y term ^vHfiijIaz^y^ . Dimensionless 2;- integration gener- 
ates a factor \meaz^uj\. Integration with respect to y adds the factor [/(i^)] 
in Eq. ffTOj) . where 



m 



1 + 



A2 i-2 



1 + 



1 {L-2)\ (hL L + 2 

'ZlI 



(11) 



is obtained by conventional methods (IGradshteyn fc Ryzhikl . Il965l ) and af- 
ter some simple algebra. Here = ^/h/eBo is the 'magnetic length' 
( lAharonov fc Bohnj . Il959l ). The L-dependence in the brackets vanishes for 
Landau levels with L < 2. The two lowest order Landau levels experience 
a constant energy shift. The final result for the Landau energy shift is the 
approximate expression given on the right in Eq. f lTOj) . Since a > (being 
the second derivative of B{s) in the minimum) is always positive, the shift in 
perpendicular energy caused by the bounce perturbation on a Landau level 
is positive and small. This follows from the rapid decay with increasing L of 
the negative term in parentheses. In general, the smallness of the correction 
is due to the pre-factor in the second term in f{L); the correction decays the 
higher the Landau level. 

To first and higher orders the bounce also introduces a weak z-dependence 
of the perturbed Landau wave function which destroys the parallel degener- 
acy. Its exact expression, taking care of the degeneracy, can be calculated by 
conventional perturbation methods, being here, however, of lesser interest as 
long as no practical application is urgent. In such for instance in 

the case of a very strong magnetic dipole field like in pulsars or magnetars), 
one would rather include the correct functional s-dependence of the magnetic 
field. Here we have just shown the expected dominant effect in the lowest 



order energy shift when electrons bounce in a magnetic mirror field geometry. 
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